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We have calculated ground state properties and excitation spectra for Ce metal with the ab 
' initio computational scheme combining local density approximation and dynamical mean-field the- 

, ory (LDA+DMFT). We considered all electronic states, i.e. correlated /-states and non-correlated 

S-, p- and d-states. The strong local correlations (Coulomb interaction) among the /-states lead 
to typical many-body resonances in the partial /-density, such as lower and upper Hubbard band. 
Additionally the well known Kondo resonance is observed. The s-, p- and d-densities show small to 
mediate renormalization effects due to hybridization. We observe different Kondo temperatures for 
a- and 7-Ce {TK.a ~ 1000 K and rA',7 ~ 30 K), due to strong volume dependence of the effective 
hybridization strength for the localized /-electrons. Finally we compare our results with a variety 
of experimental data, i.e. from photoemission spectroscopy (PES), inverse photoemission spec- 
^ , troscopy (BIS), resonant inverse photoemission spectroscopy (RIPES) and magnetic susceptibility 

I ' measurements. 
U. 

C/5 , 71. 27. -fa Strongly correlated electron systems , 74.25.Jb Electronic structure 

. 

■ Ce metal is the simplest lanthanide compound with only one atom in a face centered cubic (fee) crystal structure 
and a relatively small set of relevant electronic states derived from s-, p- d- and /-orbitals of Ce. It shows an unique 
I . isostructural (fee to fee) a 7 phase transition with increasing temperature. The high-temperature 7 phase has 15% 
larger volume and displays a Curie- Weiss-like temperature dependence of the magnetic susceptibility signaling the 
existence of local magnetic moments while the a-phase has a Pauli-like temperature independent paramagnetism [0 . 
^ ■ While many different models were proposed to describe this system (for a review see Q), the most relevant seems 
to be the periodic Anderson model. Studies based on the single impurity Anderson model Q with a hybridization 
function obtained from LDA band structure calculations were rather successful in reproducing Kondo scales and 
spectra for a- and 7-Ce. However, an empirical renormalization of the hybridization function and position of the 
impurity level were needed for satisfactory agreement between calculated and experimental spectra. 

Due to the recent development of the Dynamical Mean-Field Theory Q a more realistic treatment of Ce is now 
possible. In contrast to the Hubbard model (degenerate and non-degenerate), where hybridization occurs only between 
' correlated d- or /-orbitals, Ce is much more complicated. The direct /-/ hybridization is of the same order of 
\ magnitude as the hybridization of /-orbitals with the delocalized spd-states. Thus in order to describe Ce one 
' even has to go beyond the periodic Anderson model, where only hybridization of the correlated /-orbitals with the 
— , ] delocalized states is included. In order to address this problem we used the most general procedure for calculating 

. the Green function using a full basis set {s,p,d,f) Hamiltonian with the integration over Brillouin zone in k-space. 
^ ■ Recently a LDA-I-DMFT approach was proposed, with different methods to solve the DMFT equations: IPT 

NCA and QMC @-[l2). All these three methods were used to investigate Lai-^^Sr^rTiOa [||J|,0. The same 
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strategy was formulated by Lichtenstein and Katsnelson as one of their LDA++ approaches. Lichtenstein and 
, Katsnelson applied LDA+DMFT(IPT) Q, and were the first to use LDA+DMFT(QMC) Q to investigate the 
O ■ spectral properties of iron. Liebsch and Lichtenstein also applied LDA+DMFT(QMC) to calculate the photoemission 
spectrum of Sr2Ru04 p^ . 

Here we present results obtained within LDA+DMFT (NCA) ||,^. After convergence of the DMFT selfconsistent 
loop the hybridization function corresponds to fully interacting effective medium for the /-shell of a Ce ion. The 
interacting Green function obtained in our calculation allows us to compute ground state properties like orbital 
occupation values, Kondo temperatures Tr-, magnetic susceptibility x(0), as well as excitation spectra. Our results 
show a satisfactory agreement with the experimental data without using any adjustable parameters. 

In the following we will concentrate on a simplified local interaction. Here, we introduce two distinct Coulomb 
parameters: the intra-orbital Coulomb energy U has to be considered in case of a doubly occupied orbital, while 
the inter-orbital Coulomb energy U' applies for example in the case of a doubly occupied /-shell with electrons on 
/-orbitals with different indices. Since we neglect any exchange correlations, which is typically of the order of one 
tenth of the Coulomb interaction, we chose U = U' in order to fulfill the condition of rotational invariance of the 
local interaction UM. A more sophisticated interaction term has been already investigated within the framework of 
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LDA+DMFT(NCA) scheme We thus arrive at an interaction of the form 

HI.°';''J ^U'^nni^nrnl + Y f^rnaflm' a' ■ (1) 

The most important feature of the DMFT is that the proper one-particle selfenergy due to the local Coulomb inter- 
action is purely local Q. Thus, we obtain as an expression for the full Green function of the interacting system 

G(^) = W- ^^'^ S(z)Y^)-i , (2) 

k 

noninteracting one-particle hamiltonian h (fc) and consequently G{z) will in general be matrices in orbital space, I f 

is the diagonal matrix with matrix elements equal 1 for f-orbitals and zero for all others, I is unit matrix. The 
/c-summation is done by a standard tetrahedron method Within this method one can easily treat hybridization 
effects between correlated and non-correlated states. 

As starting point of our calculation we determined the one particle LDA Hamiltonian with the LMTO method |l^] 

considering the 6s,6p,5(i and 4/-shells. The noninteracting one-particle hamiltonian h {k) was obtained by subtracting 
the Hartree contribution of (Q) from the LDA results in order to avoid double counting The value of the 

Coulomb interaction was calculated by a supercell method and found to be C/ w 6 eV . The chemical potential was 
adjusted to conserve the number of paricles (4 electrons per site) during the selfconsistent LDA+DMFT calculation. 
Analizing the partial densities of states one can observe for both a- and 7-Ce at a temperature of T = 580 K 
intermediate size renormalization effects, in particular broadening and shifts of structures, for s- and d-states and 
only marginal effects for the p-states. This is a consequence of hybridization, which is seen by a non-vanishing s- or 
d-density at the position of the /-states in the LDA result. The /-states are strongly renormalized. A lower Hubbard 
band (LHB) is observed at the position of the corrected (double counting) LDA /-state (at about —3 eV). The 
upper Hubbard band (UHB) is situated at about 4.5 eV and describes an excitation of a doubly occupied /-state. 
At the Fermi level one observes a Kondo resonance (KR) for a-Ce, which can be described by a singlet formation 
between an unpaired /-electron and the surounding conduction electrons. The d-density for a-Ce shows an onset 
of a hybridization gap, which is well known in model calculations for the periodic Anderson model |l9|] and is a 
consequence of the formation of the singlet state between the unpaired /-spins and the conduction electrons. For 
7-Ce one observes only the onset of a KR as a consequence of the smaller compared to the a-phase and thus no 
hybridization gap opens in the conduction electron density. 

In Table | we show a comparison of our results with the results of spectral fits to electron |^ and high-energy 
neutron spectroscopy pO| ]. The occupation probabilities are in very good agreement with the experiment, as well 
as the number of occupied /-states per site, Uf. The occupation probabihties Pq, Pi and P2 for the states /", /^ 
and p were calculated from the particular ionic propagators. The Kondo temperature for a-Ce Tr-,q w 1000 K is 
roughly given by the width of the KR. Since for 7-Ce the Abrikosov-Suhl resonance is not yet well developed, it makes 
no sense to estimate Tk from the width. Instead, for 7-Ce we estimate the Tr- from the ratio of the hybridization 
strength at the Fermi level of both considered materials (Q{AQ(£i;')}/5J{A-y(£p')} = 2). With this relation we obtain 
7^,7 ~ -^Tk ci dfj- Since the NCA for multi-band models typically underestimates Tr-, it is obvious, that our absolute 
values cannot be expected to match with the experiments. Nevertheless, the ratio of the Kondo temperatures for the 
two different phases should be meaningful and are in good agreement to the experiment (see Table ^. The static 
susceptibilities x(0) calculated with the Tk as an input via x(0) = ^ , naturally have the same qualitative 

character. In this formula C is the Curie constant for the lowest 4/-state with J = | and n{p) is the occupation of 
the /° states. Again, the ratio of the susceptibilities for the both phases are in good agreement with experiment (see 
Table ^ . Also in Table || the parameter Aa„ represents the averaged value of the imaginary part of the hybridization 
function 5{A(ti;)} for an energy interval from —3 eV to eV j22| . It is also in reasonable agreement with experiment. 

For a-Ce the evolution of the KR is clearly observed down to a temperature of T = 580 K, whereas in 7-Ce the 
evolution of a quasi-particle resonance is strongly suppressed in this temperature regime by a smaller hybridization. 
Nevertheless by decreasing the temperature down to T = HQ K the onset of the many-body resonance can be observed 
for 7-Ce, too. 

A comparison of the imaginary part of the hybridization function for both phases leads to the result that a strong 
renormalization takes place in comparison to our pure LDA results, which are in agreement with Ref . . The total 
weight of this quantity for a-Ce is twice as high as the one for the 7-phase. Here it turns out that the consideration 
of all hybridization processes is extremly important in order to get a qualitatively good agreement to experimental 
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results. If one uses a model with correlated /-states only, one has to consider one electron per per site. This would 
produce a /-density with the Fermi energy in between the /^ and /^ charge excitation peak, the LHB and UHB 
respectivly. Since the hybridization function is proportional to the density of states one would observe only a small 
hybridization at the Fermi energy. Thus the additional s-, p- and d-conduction states strongly contribute to the 
hybridization function at the Fermi energy and lead therefore to different Tk- 

In the PES data for a-Ce in upper part of Fig. |^ the observed peaks are identified as /-contributions to the density 
of states by a cross section argument using different photon energies . Thus we compare the experiment with the 
calculated partial density for /-states. The theoretical /-spectrum shows a LHB which is also seen in the experiment. 

The BIS spectrum for a-Ce shows a main structure between 3 eV and 7 eV, which is attributed to Af^ final 
state multiplets. In the calculated spectrum all excitations to 4/^ states are described by the featureless UHB. As a 
consequence of the simplified interaction model all doubly occupied states are degenerate. This shortcoming in our 
calculation is responsible for the sharp peaked structure of this feature. The neglected exchange interaction would 
produce a multiplet structure, which would be closer to the experiment. The experimental peak at about 0.5 eV is 
attributed to two 4/^ final states, which are split by spin-orbit coupling. The calculated /-spectrum shows a sharp 
KR slightly above the Fermi energy, which is the result of the formation of a singlet state between /- and conduction 
states. We thus suggest that the spectral weight seen in the experiment is a result of this KR. Since we did not 
yet include spin-orbit coupling in our model, we of course cannot observe the mentioned splitting of the resonance. 
However, as it is well known 124], the introduction of such a splitting would eventually split the KR. If we used the 
experimentally determined value of about 0.3 eV for the spin-orbit splitting |25| , the observed resonance of width 
0.5 would indeed occur in the calculations. 

In the lower part of Fig. |l| a comparison between experiment and our calculation for 7-Ce is shown. The most 
striking difference between lower and upper figures is the absence of the KR in the high temperature phase (7-Ce; 
transition temperature 141 K [0) which is in agreement with our calculations. 

In Fig. 1^ our results for the non-occupied states in the /-density are compared with RIPES data |^ . The calculated 
/-spectra were multiplied by the Fermi-step function and broadened with an Lorentzian of the width 0.1 eV in order 
to mimic the experimental resolution in the theoretical curves. Here, as above the theoretical overestimation of the 
UHB is a consequence of the simplified local interaction and thus of the missing multiplet structure of the 4/^-final 
states. The main feature of the experimental spectra: strong decreasing of the intensity ratio for KR and UHB peaks 
going from a- to 7-phase, can be also seen for theoretical curves. 

In conclusion we have described a realization of a combination of density-functional theory in the local density 
approximation and the dynamical mean field theory to obtain a first-principles computational scheme for Heavy- 
Fermion systems. The scheme was set up for the first time with a combination of correlated and non-correlated 
states in order to introduce the important effect of hybridization between s-, p-, d-states and strongly correlated 
/-states. The solution of the DMFT equations was done by using the Non-Crossing approximation. We calculated 
the one-particle spectra for a- and 7-Ce and found Kondo temperature values {TK,a ~ 1000 K and Tr-^^ « 30 K), 
which explain the experimental results. 

We observe quite resonable results concerning occupation probabilities Pq, Pi, P2 and the number of /-electrons per 
site n/. The ratio of Tk and thus the static susceptibilities x(0) values for two phases are in fair agreement with the 
experimental results considering the problems of the NCA method. Moreover we found qualitative good agreement 
with PES, BIS and RIPES experiments, i.e. the position of LHB, UHB and the Kondo resonance. 

This work was partially supported by the DFG grant PR 298/5-1&2 and Russian Foundation for Basic Research 
Grant No. RFFI-98-02-17275. 
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FIG. 1. Comparison between combined PES and BIS experimental (circles) and theoretical (solid line) /-spectra 
for a- (upper part) and 7-Ce (lower part) at T = 580 K. The relative intensities of the BIS and PES portions are roughly 
for one 4/ electron. The experimental and theoretical spectra were normalized and the theoretical curve was broadened with 
resolution width of 0.4 eV. 
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TABLE I. Comparison between LDA+DMFT(NCA) calculated parameters for both a and 7 phases at T = 580 K and 
experimental values. 
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